Critical exponents in spin glasses : numerics and experiments 
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We give an overview of numerical and experimental estimates of critical exponents in Spin Glasses. 
We find that the evidence for a breakdown of universality of exponents in these systems is very 
strong. 
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I. INTRODUCTION 

It is well known that at a continuous phase transition, 
striking critical behaviour is observed. As the transition 
temperature is approached from above or below, there 
are power law singularities in a number of physical pa- 
rameters ( the specific heat, the susceptibility,...). On 
very general grounds it can be shown that various criti- 
cal exponents which govern the singularities are related 
to each other through scaling relationships. Even more 
remarkable is the fact that systems which are very dif- 
ferent from each other at the microscopic level can be 
arranged into universality classes : within a given class 
all members have strictly identical exponents. Classes are 
defined by a restricted number of parameters - basically 
the space dimension and the number of components of 
the order parameter of each system. For standard second 
order transitions this behaviour can be fully understood, 
and the exponents calculated a priori thanks to renor- 
malization group theory. This is one of the outstanding 
achievements of statistical physics (see Fisher Q| for an 
enlightening historical survey). 

At standard ferromagnetic transitions, the exponents 
follow mean field behaviour down to dimension 4 (the 
upper critical dimension). At dimensions d lower than 
4, the exponents have been calculated though the e ex- 
pansion, where e = 4 — d. There are well behaved series 
in increasing powers of e which allow one to give renor- 
malization group estimates of the exponents. The values 
can be compared with those measured from high pre- 
cision numerical simulations, and contact can be made 
with the exact values at dimension 2 (where for the Ising 
system, the exponents have been known exactly for more 
than 50 years ||). Renormalization group values, nu- 
merical values, and experimental values are all in ex- 
cellent agreement with each other. The only exponent 
whose value is not quite so thoroughly established is the 
dynamic exponent z ; even here disagreements between 
different estimates are small. Universality is a general 



rule in systems with standard second order transitions, 
except for a restricted class of two dimensional systems 
with conflicting interactions. Here it has been shown an- 
alytically as well as numerically that the exponents vary 
continuously as a function of the ratio of the competing 
interactions Hf|. This behaviour can be explained for 
these particular systems in terms of marginal operators 
in the renormalization scheme ||. 

In the Spin Glass (SG) context, for a long time it was 
by no means obvious that there were well defined phase 
transitions at all in real three dimensional materials, or 
even in 3d model systems. Although the cusp temper- 
ature is clearly marked experimentally, the specific heat 
shows no visible singularity, and the susceptibility does 
not diverge in the region of the cusp temperature. A very 
important step forward was the realization in the in the 
late 1970s that the appropriate parameter to measure is 
not the standard linear susceptibility but the non-linear 
susceptibility ra. The Miyako group in Sapporo in pio- 
neering work p| showed that there is a divergence of \ni 
at the cusp temperature ; this work and that of other 
groups J^dT^] convinced the community that there is in- 
deed a bona fide transition in a SG. With the existence of 
a transition established, estimates were found of the crit- 
ical exponents. Numerical work soon followed |]l3]-|l5|. It 
was important to measure these values in order to com- 
pare with the renormalization group approach. Does the 
combination of frustration and randomness which char- 
acterize SGs modify the basic physics of transitions in a 
fundamental way or not ? This has turned out to be a 
long story, which is still not finished. Unfortunately, as 
no clear theoretical guidelines appeared, the enthusiasm 
for the subject dropped and even the empirical ground 
rules were not fully established. Our aim here is to show 
that this study, both numerical and experimental, is well 
worth pursuing. 

Before looking in detail at the data, we can first note 
that the exponents which come out of the numerical or 
experimental analyses are very different from the clas- 
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sical ferromagnetic values. For the latter, in dimension 
3 (for either Ising, XY, or Heisenberg spins) the expo- 
nents a and r\ are numerically small, and z is near 2. For 
the spin glasses a is strongly negative (around —2 which 
is consistent with the lack of a visible singularity in the 
specific heat )and r\ is far from zero ; experimental val- 
ues range from +0.4 to —0.5 and simulation values are 
always distinctly negative, z is strong - generally around 
6. The major differences compared with the standard 
second order transition values already indicate that the 
spin glass transition lies in a quite different category. 

In Ising Spin Glasses (ISGs) the upper critical dimen- 
sion is 6 and one could imagine that from dimension 6 
down, a similar renormalization group approach would 
be valid mutatis mutandis as for the ferromagnets. In 
fact things are much more complicated. The e expansion 
has been calculated to order three |[(|, but the succes- 
sive factors in increasing powers of epsilon grow rapidly 
so it is not at all obvious how the total series will sum. 
Valiant efforts have been made over many years to set 
the theory on a firm footing using field theory ||l7| , but 
so far the only clear result is to confirm that the leading 
term in 6 — d from the e expansion is correct. After that, 
numerous badly controled terms proliferate and theory 
is of little practical help in predicting exponents even at 
d= 5. 

As a predictive theory is lacking, we are forced to turn 
to numerical and experimental methods so as to estab- 
lish the empirical values of the exponents. The empirical 
results show clear violation of the Universality rules. 

II. NUMERICAL RESULTS 

In the following discussion we will concentrate mainly 
on the systems which have been studied the most fully 
- ISGs on (hyper)cubic lattices with random unbiased 
interactions between near neighbours. 

The definitions of the critical exponents are familiar, 
with appropriate modifications for spin glasses to take 
into account the fact that that the order parameter is 
the Edwards Anderson parameter. The specific heat ex- 
ponent is a, the order parameter exponent below the 
ordering temperature is /3, the spin glass susceptibility 
exponent is 7. The non-linear magnetization at the or- 
dering temperature exponent is 5, with 

M nl ~ H 2 ' s (I) 

. The correlation length exponent is is, and the expo- 
nent for the form of the correlation function at the or- 
dering temperature is rj. The relaxation time dynamic 
exponent is z. The scaling relationships between these 
exponents are a + 2/3 + 7 = 2, dv = 2 — a, 7 = (2 — rf)v, 
v{d - 2 + if) = 2/3 and 5 = (d + 2 - rf)/{d - 2 + rf). Nu- 
merically each of the exponents can in principle be mea- 
sured independently through temperature dependences 



on large samples , though frequently they are measured 
through using finite size scaling relationships. Experi- 
mentally, 7, a, S and the combination zv can be measured 
directly while f3,v and r\ can only be obtained through 
scaling. There are many other useful relationships ; for 
instance the relaxation of the autocorrelation function at 
T g has an exponent 

x= (d-2 + r])/2z (2) 

In addition to the standard critical exponents other 
exponents can be defined, in particular the stiffness ex- 
ponent 9. This is not a critical exponent but is defined by 
the size dependence of changes in energy with boundary 
conditions [jl8|| . For a spin glass energy measurements 
can be made with periodic and anti-periodic boundary 
conditions. The variation of the sample to sample fluc- 
tuations of the energy differences scale as L e . If the zero 
temperature 9 is positive, then the ordering temperature 
is greater than zero. 

A first type of quasi-numerical approach is furnished 
by high temperature expansions |lj|,[l9],^0| . The method 
consists of an extrapolation from a finite number of exact 
terms in the high temperature series expansion of some 
thermodynamic function to its asymptotic coefficients. 
The asymptotic form of the series contains the informa- 
tion on the singularities of the function. The extrapola- 
tion is not exact, but excellent results have been obtained 
in regular systems. The situation is less favourable in 
disordered systems. Careful analysis of the spin glass 
susceptibility from a series with a large number of terms 
(up to 20 in dimension 3) provides a set of estimates 
for the values of T g and 7 obtained from different ap- 
proximant functions. If the series expansion was infinite, 
the method would become exact but in practice the lim- 
ited length of the series means that the estimates are not 
perfect. The exponent results become less accurate as 
one gets further from the upper critical dimension. For 
dimensions 5 and 4 the longest series give high quality 
estimates which can be used as independent yardsticks 
to compare with the Monte Carlo data which we will dis- 
cuss below. The method does not have the same prob- 
lems (such as thermalization) which are encountered in 
Monte Carlo simulations, but considerable know-how is 
needed to calculate long series and to extract reliable ex- 
ponent estimates from the raw series results. Up to now 
the only specific ISG series expansion results published 
are for binomial (± J) interactions. 

For dimension 5 a reliable and accurate ordering tem- 
perature and set of exponents is given by pQ | . For di- 
mensions 4 and 3 the results quoted in p9[ are more 
transparent ; it is clear that the approximant data points 
cluster satisfactorily with a strong correlation between 
estimates for T g and those for 7. 

The most widely used technique for determining criti- 
cal exponents numerically has been that of Monte Carlo 
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simulation. Many efforts have been made to measure ISG 
critical exponents accurately, despite considerable tech- 
nical difficulties. For measurements exploiting the finite 
size scaling method [O, each sample is first annealed 
numerically until the spin system can be judged to be in 
thermal equilibrium ; then the fluctuations in the auto- 
correlation function q(t) =< Si(0)Si(t) > are measured. 
The precautions necessary are described in [Q. Long 
enough times must be used for each part of the procedure, 
and the time scale defining " long enough " depends on 
the size L and the temperature T. The larger the sample 
and the lower the temperature the longer the time scales, 
so it becomes very difficult to obtain significant data on 
large samples at low temperatures. Sophisticated update 
methods have been developed which aleviate this prob- 
lem to some extent 21 1. One must be sure that mea- 
surements have been done over a sufficient number of 
independent samples. Even if numerically high quality 
data has been obtained, there may be intrinsic correc- 
tions to finite size scaling which mean that the scaling 
rules (which are valid asymptoticly for large sizes) may 
not yet hold exactly for the range of sizes studied. 

An important parameter which can be deduced from 
the equilibrium fluctuations is the spin glass suscepti- 
bility XSG =< Q 2 >j directly related to the non-linear 
susceptibility. The finite size scaling form is 



X SG = L z -y{L'/»{T~T g )) 
Precisely at T g , as a function of size 
Xsg(L) ~ L 2 ~ v 



(3) 



(4) 



T g can be determined quite accurately as the highest 
temperature at which log(xsc/ L 2 ) varies linearly with 
log(L) up to large L. 

These expressions ignore corrections to finite size scal- 
ing. There should be a further factor so that for in- 
stance the spin glass suceptibility is multiplied by [1 — 
L- w j L (L>l v {T - T g )) + ...]. The correction to scaling 
exponent w has been estimated at around 2.8 in the bi- 
nomial ISG in dimension 3 p2fl , as against 0.9 for the 3d 
Ising ferromagnet and 1.6 for 3d site percolation [p3[ . It 
turns out that the strength of the correction to scaling 
can change dramatically from one system to another. 

Many authors estimate the critical temperature T g 
through the Binder the cumulant method juj. The cu- 
mulant for a given T and L is defined by a dimensionless 
combination of moments of the autocorrelation fluctua- 
tions averaged over a large number of samples : 



g L = 1/2(3- < q 4 > I < q 2 > 2 ) 



(5) 



for Ising spins. This cumulant is defined so as to go 
from zero for a high temperature Gaussian distribution 
of q(t) values, to 1 for a unique low temperature state. 
(Other related cumulants can be defined). For a con- 
tinuous transition with a critical temperature T g , gL(T g ) 



should be independent of size L, with values fanning out 
as a function of L above and below T g . Once T g has 
been established accurately, the exponents can be esti- 
mated by plotting the whole set of data for gh{T) and 
for xSG(L,T) in an appropriate scaling form. 

In practice, sample to sample fluctuations in < q 4 > 
are strong so very large numbers of samples must be mea- 
sured (the lack of self averaging is very much worse for 
< q 4 > than for < q 2 >). The crossing point can be ill 
determined because the <?l(T) curves do not fan out ap- 
preciably at temperatures lower than the crossing point. 
When this is the case, relatively minor corrections to fi- 
nite size scaling can modify the apparent crossing tem- 
perature. (Correction factors of the form given above 
should apply to both < q 4 > and < q 2 >). Finally, the 
values of the exponents deduced from the scaling plots, 
strongly correlated with the T g estimate, frequently vary 
very steeply with the apparent value of T g . In conclusion, 
the exponent estimates obtained in this traditional way 
must be treated with considerable caution, even when 
large scale numerical efforts have been made. Accurate 
and reliable results can be obtained only in favourable 
circumstances. 

When calculations have been made to large sizes, it is 
possible to use the scaling rules for the spin glass suscep- 
tibility and for the spin correlation length to estimate the 
critical temperature and the exponents by extrapolations 
to infinite size [2^ ]. This method should be reliable, as 
corrections to finite size scaling are kept well under con- 
trol. 

Alternative techniques which have been less widely 
used rely (at least partly) on dynamic measurements. In 
massive simulations on large samples (64 3 spins) Ogielski 
p"5[ studied the relaxation of q(t) as a function of tem- 
perature. An advantage of this method is that a strict 
thermal equilibrium state is not necessary ; as long as the 
anneal has been made over a time r a much longer than 
the subsequent measuring time over which q(t) is studied, 
the measured q(t) curve will be the true thermal equilib- 
rium form (see for instance ^5|). Ogielski assumed the 
standard critical behaviour for q(t), which is 



q{t) = t-*f(t/r(T)) 



(6) 



He estimated T g from the divergence of the relaxation 
time t(T). With T g in hand he estimated the critical ex- 
ponents from a combination of dynamic and equilibrium 
measurements. 

One can also exploit the critical behaviour of strictly 
non-equilibrium dynamics. Suppose a spin glass sam- 
ple is initially at infinite temperature (so the spins have 
random orientations) ; it is then annealed from this 
configuration to the critical temperature T g . The non- 
equilibrium spin glass susceptibility at a time t after the 
start of the anneal will increase as t h fE6|where 



h = (2- v )/z 



(7) 
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Analagous non-equilibrium dynamic parameters have 
been studied extensively in a large number of regular 
systems pTj . This non-equilibrium scaling behaviour has 
been established on a very firm theoretical basis (2£| . An 
obvious practical advantage for numerical work is that 
no preliminary anneal is required. 

Now it is clearly possible to combine the measurements 
of the dynamic relaxation exponent x(T) and h(T) at a 
series of test temperatures Ti to obtain a sequence of 
apparent exponents rj(Ti) and z(Ti) : 



For the Gaussian distribution, from the values of 
the stiffness exponent 9(d) as a function of dimension 
Jl8|,^l 33 1, we can estimate the lower critical dimension 
di , the dimension at which 9 = 0, figure [l| Because the 
Gaussian distribution is continuous, as di is close to 2.50, 
ij(di) must be close to —0.50, from equation (10). 



r)= (4x + (2-d)h)/(2x + h) 



(8) 



z = d/(2x + h) 



(9) 



The rj{Ti) can be compared with independent rj(Ti) es- 
timates from the equilibrium spin glass susceptibility 
(equation (9)). Consistency dictates that at the true T g 
the two estimates must coincide. This leads to estimates 
of T g , 77 and z which turn out in practice to be precise 
(the crossing is clean ) and virtually free from pollution 
by corrections to finite size scaling [^9). Once these pa- 
rameters are well established, the other exponents such 
as v can be determined from conventional scaling plots 
with only one unknown parameter. We will refer to this 
method as the " three scaling rule " technique. 

The low dimensions present special cases. First, in 
dimension 1 where critical temperatures are certainly al- 
ways zero for short range interactions, the stiffness ex- 
ponent 9 is exactly —1 for continous distributions of in- 
teractions and for ±J interactions |l8|. In dimension 
2 it is also well established that the ordering tempera- 
ture T g = |Q (except possibly when the interactions 
are ±J ]3(|). From the definition of 77, when T g is zero 
and for a unique ground state (corresponding to any con- 
tinuous distribution of interactions such as the Gaussian 
distribution) there is an additional scaling rule ||l4| 



77 = 2 - d 



(10) 



As fewer spins are involved for a given size L than at 
high dimension, it is easier to cover a wide range of L 
for finite size scaling in Monte Carlo simulations. Even 
better, sophisticated numerical techniques exist to find 
exact ground states for systems up to large sizes [ j33| . 
(Recently exact ground states in dimension 3 have also 
been obtained to quite large L fl32|). 



III. EXPONENT VALUES 

We will concentrate on the values of 77, and start from 
the low dimensions. 

For the 2d binomial ISG, 77 has been estimated to be 
0.20 ± 0.02 Q. Curiously, even if the T g is small but 
non-zero, the estimated value of 77 is very similar p0|. 



FIG. 1. Stiffness exponent 9 as function of dimension for 
ISGs with Gaussain interactions. Dimension 1 j^], dimension 
2 dimension 3 H) . 

Numerous estimates have been given of the exponents 
for ISGs in dimension 3 with binomial or Gaussian inter- 
action distributions. We have summarized the situation 
in p2[ , where we find that there are strong deviations 
from finite size scaling for the binomial case and where 
we obtain reliable exponent estimates principally from 
the three scaling rule method. Recent data on the bino- 
mial case established by extrapolation to infinite size |24|] 
are consistent with the values of Other estimates 

relying on the Binder cumulant method are probably af- 
fected by corrections to finite size scaling and so are less 
reliable. 

We will use the results obtained in dimension 4 to 
demonstrate the coherence of the different methods when 
these are used carefully, and the inescapable conclusion 
that the values of the exponents change with the form of 
the interaction distribution. 
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FIG. 2. Binder cumulant data for the binomial (± J) inter- 
action ISG in dimension 4, (A. P. Young, unpublished). Sizes 
L — 4, 6, 8, 12. The curves should intersect at the critical 
temperature. 

Figure || shows high precision data for the Binder cu- 
mulant in the binomial case p4| ]. There is a clean inter- 
section of the curves at T = 2.00 ± 0.01. 
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FIG. 3. r/(T) data for the binomial interaction ISG in di- 
mension 4. r/i (squares) is from the spin glass susceptibility 
(A. P. Young, unpublished) and r\2 (circles) is from the h and 
x effective exponent, see text. The intersection should corre- 
spond to the critical temperature T g and r)(T g ). 

Figure || shows the intersection of the two curves for 
n{Tj) used in the three scaling method described above 
|35| . The intersection is at precisely the same temper- 
ature, validating this technique. From the intersetion 
point we can deduce T g and 7/ = 0.30 ± 0.02. As the 
slope of the curve for 77 (T^) for the estimate from the 
dynamic exponents x and h is much weaker than from 
that deduced from the equilibrium finite size scaling spin 
glass suceptibility, the precision on the value of 77 is much 
higher using this method rather than working only with 
the pure equilibrium susceptibility and Binder cumulant 
data. The series expansion data |l9|] is in quite indepen- 
dent agreement with these Monte Carlo results. ^From 
the results plotted in Figure 3 of for this value of 
T g one would expect 7 = 2.1 ± 0.1, or v = 0.92 ± 0.05 
from the scaling relation. This value of v is perfectly con- 
sistent with direct scaling of the Monte Carlo spin glass 
susceptibility data. 
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FIG. 4. Binder cumulant data for the Gaussian interaction 
ISG in dimension 4. Sizes L — 4,6, 10. 

Figure |I] shows Binder cumulmant data for the 4d 
Gaussian case ; it can be seen that the intersection point 
lies at T — 1.785 ±0.01. Figure || shows the three scaling 
rule intersection. Again the agreement is excellent, but 
the value of rj at the intersection point r\ = —0.44 ± 0.02, 
is considerably higher than for the binomial case. 
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FIG. 5. ri(T) data for the Gaussian interaction ISG in di- 
mension 4. Squares correspond to the SG susceptibility, cicles 
to the (h, x) method. The r\ value at the intersection is clearly 
different from that of the binomial interaction distribution. 

Figure || shows a direct plot oflog(xsc/ L 2 ) against 
log(L). This has to be straight at T g . It can be seen 
that T g must lie just below T = 1.79, and rj can be esti- 
mated from the corresponding slope. We can note that 
there is excellent agreement point by point between the 
data of pf^ ] and the present results wherever comparisons 
can be made. Using a scaling plot we find v = 1.08±0.10. 
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FIG. 6. The spin glass susceptibility \sg normalized by 
L 2 against L on a log-log plot for the Gaussian interaction 
ISG in dimension 4. Temperatures from top to bottom are 
1.75,1.77,1-79, 1.811.83,1.85. 

The 4d binomial exponent values given here are consis- 
tent with but more accurate than those obtained in jj7) . 
For the Gaussian system, the value of r] = —0.35 ± 0.05 
quoted in Q is low because of a marginally overesti- 
mated Tg, 

A careful analysis of the data on these two systems in 
dimension 4 but with different interaction distributions 
shows conclusively that the exponents are different, so 
universality does not hold. 

In dimension 5, there is excellent agreement between 
series estimates |2(| and Monte Carlo estimates [^5) for 
the binomial system. 



upper critical dimension value r] = 0. The straight line 
starting from this point is the leading term in the e ex- 
pansion : T](d) = —(6 — d)/3 + ... |l6|. It is clear that 
the data demonstrate that the exponents for the two sys- 
tems, binomial and Gaussian, follow regular curves as a 
function of dimension. Independent results are consistent 
with each other. However the values at each dimension 
are different for the two sets of interactions so universality 
is clearly violated. Results for other sets of interactions 
p9| , p5[ confirm this variation. 

Critical exponents have also been measured in the Ising 
phase diagram as a function of the ratio p of the number 
of ferromagnetic to antiferromagnetic bonds. In the fer- 
romagnetic region to the right of the Nishimori line, the 
static exponents do not appear to vary with p, but the 
dynamic exponent z does change continuously and very 
significantly |3S|| . 

We can note that the Migdal-Kadanoff renormalization 
approach (which should be exact for a heirachical lattice) 
has been used to measure effective ordering temperatures 
and exponents for four different ISG interaction distri- 
butions. ]3S[ |. For dimension 3, a diamond hierachical 
lattice, and a renormalization factor b — 2, the order- 
ing temperatures are in excellent agreement with Monte 
Carlo estimates on cubic lattices. However the Migdal- 
Kadanoff calculations lead to a universal saddle point 
critical temperature and exponent values. This result 
seems to be closely related to the hypotheses intrinsic to 
the Migdal Kadanoff method. 
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FIG. 7. The critical exponent r\ as a function of dimension 
for ISGs with binomial (squares) and Gaussian (circles) inter- 
actions. The upper straight line represents the leading term 
in the e expansion [[l6| starting from the upper critical dimen- 
sion 6. The lower straight line is equation (10). The points 
in dimensions 3,4 and 5 are from the "three scaling rules" 
method. The dimension 2 point is from The lower criti- 
cal dimension point for the Gaussian interactions is explained 
in the text. 

In figure |7] we display results for r\ as a function of di- 
mension for ± J and Gaussian interactions, including the 
point at di for the Gaussian case and the point at d = 2 
for the ± J case. The point at d = 6 corresponds to the 



IV. SPIN GLASS EXPERIMENTS 

Estimates of exponents have been made by many ex- 
perimental groups, in general using slightly different pro- 
tocols and on different materials. The non-linear suscep- 
tibility is defined as 

Xm = Xo - M/H (11) 

where \o is the zero field magnetic susceptibility. At T g 
Xni should behave as H 2 I S and above T g 

M = XoH + X2H 3 +XiH 5 + ... (12) 

with \2 diverging as t~ 7 , where t is the reduced temper- 
ature. The Suzuki equation of state is Q 

Xni = t^g(H 2 /t^) (13) 

Both d.c. and a.c. magnetization techniques have been 
used . For instance Monod and Bouchiat and Bouchiat 
used d.c. entirely, taking care to stay in a field range 
where the non-linear magnetization remained less than 
10% of the linear magnetization so as not to corrupt the 
data. Svedlindh et al |ll| first analysed low field and low 
frequency a.c. measurements to fit to 7 and T g . (The 



G 



ordering temperature was in agreement with the low d.c. 
field cusp temperature). They then used the equation 
of state with d.c. measurements up to moderate fields 
with (3 as the only fit parameter. In a second set of 
experiments they measured the field and frequency 
dependence of the a.c. susceptibility and obtained an 
estimate of the dynamics exponent zv. In a sophisticated 
experiment, Levy measured the Fourier transform 
spectrum of the magnetization response to a 0.1 Hz field, 
picking up a series of non-linear susceptibilities from the 
different harmonics. He could deduce accurate values of 
static and dynamic exponents. 

The experiments have to be performed with care, per- 
ticular attention being paid to the proper identification of 
the critical temperature and to the necessity to remain 
in a suitably low field range throughout. Recent mea- 
surements in which the exponent 8 was measured 
in a range of different materials using one single protocol 
gave excellent confirmation of earlier data by other teams 
(except for the case of AuFe where an early experiment 
had given values of exponents out of line with all other 
results). These results validate the earlier measurements 
and show that the considerable spread of values of ex- 
ponents reported for different materials, is not due to 
artefacts in the different measuring procedures. As in 
the numerical data it is evident that the expected uni- 
versality of exponents breaks down. 

The finite temperature T g values for real material spin 
glasses have been somewhat of a mystery for some years. 
These systems are Heisenberg, and reliable numerical 
work has demonstrated that Heisenberg spin glasses in di- 
mension 3 should have zero temperature ordering |45| , [46[ . 
Kawamura |46| has made the interesting suggestion that 
the ordering process in real life Heisenberg materials is 
basically a chiral spin glass ordering. This ordering would 
not be visible directly to magnetization experiments if 
there were no anisotropy. However in all real systems ran- 
dom anisotropy (of the DM type ^7j) is always present, 
and by coupling the chiral degrees of order with the mag- 
netism an anisotropy, however weak, reveals the chiral 
order. 

The critical exponents for pure Heisenberg chiral or- 
dering in dimension 3 have been estimated numerically 
. The best values are around v = 1.25 and r\ = +0.7. 
It can be noted that v is similar to the Monte Carlo 
Ising values which we have presented, while r] is strongly 
positive rather than negative as seen in the numerical 
Ising work. A plausible hypothesis is one in which the 
exponents change progressively from chiral-like for weak 
anisotropy to Ising- like for strong anisotropy. On this sce- 
nario, the value of v should remain relatively stable for all 
the materials, while the value of r\ should vary progres- 
sively, becoming gradually more negative as anisotropy 
increases. 

So far the qualitative trend for systems where both 
anisotropy and exponents have been measured is in ex- 



cellent agreement with this picture. All the experimental 
exponent values fall within the chiral limit at one end 
and the Ising limit at the other end. For the three alloy 
systems AgMn, CuMn and AuFe the anisotropics are 
weak, moderate and strong respectively. The v values 
are similar, near 1.3, while the r\ values are about 0.4, 0.1 
and —0.1 jlij. The trend of exponent values is clearly in 
the sense predicted by the scenario. 

V. CONCLUSION 

The main lesson which can be drawn from this 
overview of numerical and experimental exponent data 
in spin glasses is that transitions in these glassy systems 
are quite different from those in regular systems with 
standard second order transitions. The values of the ex- 
ponents are far from those in regular systems and the 
breakdown of Universality is manifest in carefully anal- 
ysed data. 

The statistical physics community has been very loth 
to accept evidence against universality because the renor- 
malization scenario appears to give such an appealingly 
general picture of behaviour at transitions. However the 
fact that it has proved extremely difficult on the field 
theory level to produce predictions for spin glasses is a 
strong indication that unexpected behaviour cannot be 
excluded a priori. 

What possible mechanism could lead to this break- 
down ? It has been found analytically that Ising spin 
systems with ferromagnetic interactions on hierachical 
lattices show no universality [[l8| . For the spherical model 
on graphs of non-integer dimension, the exponents vary 
continuously with the spectral dimension of the graph 
We can speculate that in spin glasses the effective 
dimension of the system at criticality could depend on 
the form of the interaction distribution. 

It would be unfortunate if this phenomenon was left 
unexplored because of preconceptions as to the physi- 
cal laws which should hold for complex systems. If it 
could be accepted that critical behaviour is much richer 
in glassy systems than at conventional second order tran- 
sitions without frustration, an important new field of in- 
vestigation should open out. What control parameters 
affect the exponents and why ? What are the implica- 
tions for the physics of the glass transition in the most 
general sense ? 
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